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We investigate the locality of interactions in hydrodynamic turbulence using data from a direct 
numerical simulation on a grid of 1024^ points; the flow is forced with the Taylor-Green vortex. 
An inertial range for the energy is obtained in which the flux is constant and the spectrum follows 
an approximate Kolmogorov law. Nonlinear triadic interactions are dominated by their non-local 
components, involving widely separated scales. The resulting nonlinear transfer itself is local at 
each scale but the step in the energy cascade is independent of that scale and directly related to 
the integral scale of the flow. Interactions with large scales represent 20% of the total energy flux. 
Possible explanations for the deviation from self-similar models, the link between these findings and 
intermittency, and their consequences for modeling of turbulent flows are briefly discussed. 

PACS numbers: 47.27.Eq,47.27.Ak,47.65.+a 



Flows in nature are often in a turbulent state driven 
by large scale forcing {e.g. novae explosions in the in- 
terstellar medium) or by instabilities {e.g. convection in 
the sun). Such flows involve a huge number of coupled 
modes leading to great complexity both in their tempo- 
ral dynamics and in the physical structures that emerge. 
Many scales are excited, for example from the planetary 
scale to the kilometer for convective clouds in the atmo- 
sphere, and much smaller scales when considering micro- 
processes such as droplet formation. The question then 
arises concerning the nature of the interactions between 
such scales: are they predominantly local, involving only 
eddies of similar size, or are they as well non-local? It is 
usually assumed that the dominant mode of interaction 
is the former, and this hypothesis is classically viewed as 
underlying the Kolmogorov phenomenology that leads to 
the prediction of a E{k) ^ k^^/^ energy spectrum; such 
a spectrum has been observed in a variety of contexts al- 
though there may be small corrections to this power-law 
due to the presence in the small scales of strong localized 
structures, such as vortex filaments Q. 

Several studies have been devoted to assess the de- 
gree of locality of nonlinear interactions, either through 
modeling of turbulent flows, as is the case with rapid 
distortion theory (RDT) j^j or Large Eddy Simulations 
(LES) [3 , or through the analysis of direct numerical sim- 
ulations (DNS) of the Navicr-Stokes equations (see eAi. 
0, 0, 0] ) , and more recently through rigorous bounds [a| . 
The spatial resolution in the numerical investigations was 
moderate, without a clearly defined inertial range and 
the differentiation between local and non-local interac- 
tions was somewhat limited. Thus, a renewed analysis at 
substantially higher Reynolds numbers in the absence of 
any modeling is in order; we address this issue by ana- 
lyzing data stemming from a newly performed DNS on a 
grid of 1024'^ points using periodic boundary conditions. 

The governing Navier-Stokes equation for an incom- 
pressible velocity field v, with V the pressure, F a forcing 



term and = 3 x 10 ^ the viscosity, reads: 



dv „ 



-VP + vV^^r 



(1) 



together with V • v = 0. Specifically, we consider the 
swirling flow resulting from the Taylor-Green vortex 0: 



' TG 



(fco) = 2F 



sin(A:o x) cos(fco y) cos(A:o z) 
— cos(A:o x) sin(/co y) cos(fco z) 




, (2) 



with fco = 2. This forcing generates cells that have locally 
differential rotation and hclicity, although its net helicity 
is zero. The resulting flow models the fluid between two 
counter-rotating cylinders and has been used widely to 
study turbulence, including studies in the context of the 
eneration of magnetic fields through dynamo instability 
The Reynolds number based on the integral scale L = 
2tt J E{k)k-^dk/E w 1.2 (where E is the total energy), 
is i?e = UL/v K, 4000, where U is the r.m.s velocity. 
The Reynolds number based on the Taylor scale A = 
2tt{E/ J k^E{k)dky/^ w 0.24, is Rx w 800. 

The code uses a dealiased pseudo-spectral method, 
with maximum wavenumber k^ax — 341 and kmaxV = 
1.15, where 27r77 = 27r(z/'^/e)^/^ is the dissipation scale 
and e is the energy injection rate: the flow is suffi- 
ciently resolved since 1 /rj is within the boundaries of the 
wavenumbers handled explicitly in the computation. 

Details of the flow dynamics will be reported elsewhere; 
suffice it to say that the flow reproduces classical features 
of isotropic turbulence 0: the energy spectrum is well- 
developed (see Fig. with a constant energy flux for 
k € [5, 20] and maximally helical vortex tubes are found, 
as predicted in 0| and shown in [ill IT^ . Finally, the 
anomalous exponents of longitudinal structure functions 
are in excellent agreement with previous studies up 
to order p = 8 (see Table I), including analysis without 
using the extended self-similarity (ESS) hypothesis 0|. 

To investigate the interactions between different scales 
we split the velocity field into spherical shells in Fourier 
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FIG. 1: Compensated energy spectrum and (inset) absolute 
value of the energy flux n{k) in the stationary regime. 



TABLE I: Order p and anomalous exponents (^p computed on 
two snapshots of the velocity field using the interval of scales 
with constant energy flux; the anomalous exponents Cp^^ ^.re 
computed using the ESS hypothesis. 



p 


1 


2 


3 


4 


5 


6 


7 


8 


Cp 
^p 


0.366 
0.364 


0.704 
0.699 


1.005 
1 


1.271 
1.264 


1.502 
1.495 


1.703 
1.695 


1.878 
1.869 


2.029 
2.020 



space of unit width, i.e. v = J2k^k where vk is the 
filtered velocity field with K < \k\ < K + 1 (from now on 
called shell K) 0|. Usually, octave bands are used to 
define the shells {i.e. shells of width AisT" are used, where 
AK is a constant often set to 2). This choice is based on 
the usual hypothesis that interactions are mostly local 
and self-similar in Fourier space, i.e. the nonlinear term 
in eq. couples triads of modes (fc, p, q) in Fourier space 
with k ^ p ^ q. Since we want to verify these hypotheses, 
we choose to use a linear step for the shells. This election 
does not imply any loss of generality, and if interactions 
are indeed local our results should be compatible with 
results using octave bands. 

From equation JQ), the rate of energy transfer 




K-Q 

FIG. 2: Normalized energy transfer from the shell Q to the 
shell K with Q G [10, 80]. The width of the lobes is indepen- 
dent of K and all the peaks are at — Q ~ fco . 



T3{K, P, Q) (a third-order correlator) from energy in shell 
Q to energy in shell K due to the interaction with the 
velocity field in shell P is defined as usual [islllij as: 

TsiK, P,Q) = - J VK ■ (vp • V)vQdx3 . (3) 

If we sum over the middle wave number P we obtain the 
total energy transfer T2{K, Q) from shell Q to shell K: 

T^iK, Q)^Y1 ^3(^"' P,Q) = - f^K-{v- V)vQdx3 . 
p •' 

(4) 

Positive transfer implies that energy is transfered from 
shell Q to iC, and negative from K to Q; thus, both 
T3 and T2 are antisymmetric in their (if, Q) arguments. 
T2{K,Q) gives information on the shell-to-shell energy 
transfer between K and Q, but not about the local- 
ity or non-locality of the triadic interactions themselves. 
The energy flux plotted in Fig. ^s reobtained from 
these transfer functions as n(fc) = — X]a:=o ~ 
— J2k=o X]q T2{K, Q). Note that the transfer terms de- 
fined in Eqs. (|3I4|I are integrated over all volume in real 
space. Since in periodic boundary conditions there is no 
net flux of energy through the walls, this is enough to en- 
sure that contributions due to advection (which do not 
lead to cascade of energy to smaller scales) cancel out 
(see e.g. 0). 

Figure 121 shows the energy transfer T2{K, Q) plotted as 
a function oi K — Q for 70 different values of Q varying 
from 10 to 80. For each value of Q, the x axis shows the 
different K shells giving or receiving energy from that 
shell Q. All curves collapse to a single one: the energy in 
shell K is received locally from shells with wavenumber 
K — A A' and deposited mostly in the vicinity of if -I- Ak, 
with Ax ~ ko for all values in the inertial range. In 
other words, the integral scale of the flow, related to the 
forcing scale kg^ , plays a determinant role in the process 
of energy transfer. As a result, the transfer is not self- 
similar, and the integral length scale is remembered even 
deep inside the constant-flux inertial range. 

This break down of self-similarity indicates that dom- 
inant triadic interactions can be non-local. To examine 
further this point, we need to investigate individual tri- 
adic interactions between Fourier shells by considering 
the tensorial transfer T^{K,P^Q). We will study three 
values of Q, {Q = 10, 20, and 40); for each case, P will 
run from 1 to 80, and K from Q - 12 to Q -|- 12. 

In figure 01 we show contour levels of the transfer 
T3{K, P,Q) for Q = 40. This figure represents energy 
going from a shell Q to a shell K through interactions 
with modes in the shell P. As in Fig. 12 positive trans- 
fer means the shell K receives energy from the shell Q, 
while negative transfer implies the shell K gives energy 
to Q. The strongest interactions occur with P ko, and 
therefore the large scale fiow is involved in most of the 
T2 transfer of energy from small scales to smaller scales. 
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FIG. 3: Contour levels of the transfer function T^{K,P,Q) 
for Q — 40. Solid lines correspond to positive transfer, and 
dotted lines to negative transfer. 
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FIG. 4: Comparison of the transfer functions T2{K, Q) (solid 
line), Ti'°''{K, Q) (dotted line), and TsiK, P = 3, Q) (dashed 
line), for three values of Q. 

Note that the individual triadic interactions with P ^ ko 
and JsT ~ Q ± fco are two orders of magnitude larger than 
local triadic interactions. 

When T3 {K, P, Q) in Fig. 21 is summed over all val- 
ues of P, the transfer function T2{K,Q) is recovered. 
This allow us to define the transfer rate due to inter- 
actions with the large scale flow, and due to local in- 
teractions, summing P over different ranges. Indeed, 
to further illustrate the dominance of the large scale 
flow in the involved interactions, we compare in Fig. 
21 the total transfer function T2{K,Q) with the trans- 
fer due to the large scale flow T3{K,P = 3,Q), and 
with the transfer due to local interactions in octave 
bands T^"%K,Q) ^ Ep=q/2 ^3(if, P, Q)- The figure 
indicates that the transfer due to the local interactions 
((5/2 < P < 2Q) is smaller than the transfer due to the 
integral length scale velocity field, and this behavior ap- 
pears to be stronger as the value of Q is increased. The 
remaining transfer comes from interactions with P-shells 
with wavenumbers between 1 and Q/2 (excluding P = 3), 
which are also non-local in nature. Therefore, as K and Q 
get larger (as we go further down in the inertial range), 
the dominant triads {K, P, Q) become more and more 



elongated, corresponding to more nonlocal interactions. 
As a result, detailed interactions between triads of modes 
are nonlocal, while the transfer of energy T2{K, Q) takes 
place between neighboring shells: local energy transfer 
occurs through non-local interactions. These results sup- 
port previous claims at smaller resolution 0, 0, 0] that 
a significant role in the cascade of energy in the iner- 
tial range is played by the large scale components of the 
velocity field. 

However, when computing the energy flux through a 
shell fc, i.e. integrating T2{K, Q) over all values of Q, and 
K from to /c, these non-local interactions give ~ 20% of 
the total flux, since many more local triads contribute in 
the global summation. We note that this fraction (20%) 
is independent of fc, provided that k is large enough and 
in the inertial range. 

We are left therefore with two puzzles. First, why is 
the large scale flow more effective (at the level of individ- 
ual triadic interactions) in "destroying" small size eddies 
than similar size eddies, when phenomenological argu- 
ments in the Kolmogorov spirit suggest otherwise? And 
secondly, why is the energy spectrum so close to k~^/^ in 
the constant flux region, when just advection by the large 
scale flow would suggest a shallower spectrum ~ k~^l 
(see e.g. 0]). In what follows, we give a brief review of 
possible answers as well as a simple model that shows 
how a k^^/^ energy spectrum can be obtained by advec- 
tion and stretching of the small scales just by the large 
scale flow. 

A possible answer to explain the strong non-local tri- 
adic interactions is that the Reynolds number in the 
present simulation is not high enough to observe dom- 
inance of local triads, and the decrease in amplitude of 
the small scale fields due to viscosity makes this interac- 
tions (when compared to the large scale flow) smaller. 

Another possible answer would be that the wavenum- 
ber bands defining the local interactions {i.e. the range 
of values in P used to define P;f'°'^), that were arbitrarily 
taken here to have a width of 2", could be as wide as 10" 
as some authors suggest 0. If this is the case, a DNS 
with an inertial range that spans at least three orders of 
magnitude in wavenumbers would be required to actually 
observe strong local interactions. 

However, neither of these answers address the second 
question concerning why a Kolmogorov energy spectrum 
is observed at moderate values of the Reynolds number. 
If we look at phenomenological scaling arguments, we 
see that there is one major assumption that may not be 
satisfied. Current models assume that the energy is dis- 
tributed in a hierarchy of vortices of size L, i/a, L/a^^ ... 
(with a > 1), with no specific geometry. However, ex- 
periments as well as numerical simulations have shown 
that enstrophy is distributed in vortex tubes, where two 
distinct length scales can be identified: one is the width 
of the tube / that is typically small and varies, and one is 
its length L, typically of the order of the integral scale. 
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It is not clear therefore when two such structures inter- 
act, which length scale is responsible for determining the 
time scale of the cascade. 

Prom the analysis presented here, a simple model for 
turbulent flows consistent with several features observed 
in simulations and experiments can emerge (see below). 
First, recall that Ref. found that helical vortex tubes 
capture 99% of the energy, give a k~^/^ spectrum, and 
are responsible for the strong wings in the PDF of ve- 
locity gradients. Furthermore, it was shown in Q that, 
when decomposing the velocity field in a large scale com- 
ponent U and a small scale one u, artificially dropping 
local interactions in a simulation (an operation akin to 
RDT) gives enhanced intermittency (in the sense that a 
stronger departure from linear scaling of anomalous ex- 
ponents is observed), while when non-local interactions 
are dropped the intermittency of the flow decreases ■ 

The data analyzed in the present paper implies that, at 
low order of correlators, i.e. when considering the energy 
flux, the interactions are mostly local. But when going to 
third-order individual triadic interactions (such as with 
T3), the non-local components are dominant and involve 
the integral scale. We note that this is consistent with the 
fact that departures from a linear scaling by anomalous 
exponents with the order of structure function is stronger 
as the order is increased, since it involves more non-local 
interactions linked to the geometrical structure of vortex 
tubes. This leads to a model of small-scale interactions 
involving three small scales that are substantially weak- 
ened and gaussian, thus in agreement with the findings 
in that such uw-likc terms weaken intermittency as 
well when included in the full dynamics. 

As a result, if we take into account the vortex tube 
structure of a turbulent flow, the picture of the classi- 
cal Richardson cascade may change: a possible model 
to explain the aforementioned results is to take the time 
scale of the cascade as given by the geometric average of 
the length scales involved, based on the cubic root of the 
volume of the vortex tube. If this is the case, the energy 
dissipation rate of vortex tubes with velocity ui due to the 
large scale flow Ul is given by e • Ul/ {PL)'^/'^ . This 
implies that, for constant flux, ui ~ l^/^ \J eL^I^ /Ul] this 
scaling recovers the Kolmogorov spectrum, although in a 
different spirit • Note that the spirit of this derivation 
is close to multifractal models used to explain intermit- 
tency corrections 0. 

Finally, we would like to point out that the success of 
models involving as an essential agent of nonlinear trans- 
fer the distortion of turbulent eddies by a large-scale flow 
- as in RDT and its variants or as in the alpha model 
[1^ where the flow is interacting with a smooth veloc- 
ity field (see also - may be in part explained by 
the present findings that confirm that nonlinear triadic 
interactions are mostly nonlocal and involve the integral 
scale. Similar results have alre ady been obtained for flows 
coupled to a magnetic field [l6j, where the weakening 



of nonlinear interactions may occur in different fashions, 
e.g. Alfvenization or force-free fields, and where the sec- 
ond order transfer T2 between velocity and magnetic field 
in the induction equation is itself non-local. 
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